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P3 ! We calculate the overlap between two many-body wave functions for a superfluid 

film containing a vortex at shifted positions. Comparing the results to phenomenological 

^ ' theories, which treat vortices as point particles, we find that the results are consistent 

if the point-particle vortices are considered as under the action of the Magnus force and 
> 
•^ I in weak interaction with sound waves of the superfiuid. We are then able to resolve the 

^ ■ disagreement concerning the effective mass of vortices, showing it is finite. 
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Vortices play an important role in the understanding of both static and dynamical 
properties of a superfluid[l]. They determine the Kosterlitz-Thouless phase transition [2], 
and provide a mechanism for the mutual friction between the superfluid and the normal 
fluid[3]. Due to advances in experimental techniques, there are many studies of problems 
related to vortex dynamics, such as the quantum nucleation of vortex rings induced by 
moving ions[4] and quantum phase slippage near a submicron orifice[5]. In two dimensions, 
the theoretical framework for understanding these dynamical phenomena is based on an 
effective point-particle formulation of vortex dynamics, and has been very successful [1]. 
Naturally, physical quantities in the phenomenological theory, such as the vortex mass, 
the Magnus force, and the friction should be derived from a microscopic theory. However, 
the current understanding of these quantities is in a confused state: there is no clear 
calculation of the coupling of the vortex to the low lying excitations responsible for the 
friction, and the theoretical estimates of the vortex mass range from zero[6], to finite[7], 
and to infinite[8]. There is also a suspicion that an effective mass may not be meaningfully 
defined for a vortex after all [9]. 

The purpose of the present paper is to present a conceptually straightforward calcu- 
lation to give clear constraints on these quantities. We invoke a microscopic description of 
the vortex by writing a Feynman many-body wave function for a superfluid film containing 
a vortex[10]. We calculate the overlap integral between such a state and that with the vor- 
tex shifted a distance away, and find how it behaves as a function of the distance. We also 
calculate the same quantity within the phenomenological point-vortex theory. Comparing 
the two, we conclude that the effective mass of the vortex cannot be infinite, and that the 
coupling of a vortex with low lying excitations must be sufficiently weak. At the end of the 
paper, we will discuss the generality of our approach and its application to other systems. 

Let us start with the phenomenological theory of vortex dynamics in a two dimensional 
superffuid film. A vortex is regarded as a point particle moving under the infiuence of the 
Magnus force hp^z x v, where h is the Planck constant, po is the 2-d superfiuid number 
density, z is the unit vector normal to the film, and v is the velocity of the vortex. Its 



effective Hamiltonian may be written as 

H,^^J~inV-qA{r)f, (1) 

where niy is an effective inertial mass of the vortex, the vector potential A in the symmetric 
gauge for the Magnus force is (— y, x)hpo/2, r is the vortex coordinate, and g = ±1 is the 
vorticity of the vortex. Eq.(l) can be understood by drawing an analogy with the case of 
a two dimensional electron moving under the influence of the Lorentz force by a magnetic 
held, with q interpreted as the vortex 'charge'. [11] 

This simple phenomenology is unfortunately not adequate if one wishes to compare 
with a more microscopic theory. One must also include interactions with low lying excita- 
tions such as various sound waves of the superfluid, which may be realized by the following 
model interacting Hamiltonian 

H, = qJ2 M(fc)e^''-'-(ak + al J, (2) 

k 

where k is the wave vector of a low lying excitation with the corresponding creation 
(annihilation) operator aj^ (ak)- An index labeling different kinds of excitations is omitted 
for notational simplicity. Coupling of this form conserves the total momentum of the 
system, as is necessary for a translationally invariant system. The Hamiltonian for the low 
lying excitations is 

^e = ^ hukttlai^. (3) 

k 

Therefore the total Hamiltonian of the system, a vortex and the low lying excitations, is 

H = H, + H, + H,. (4) 

In the following, we will focus our attention on the overlap integral between different 
vortex states, in which the vortex mass m„, the Magnus force, and the coupling between 
the vortex and the low lying excitations should be involved. 

In the absence of the coupling to the low lying excitations, the overlap integral between 
two coherent states for a vortex centered at Tq and tq in the ground state can be calculated 
as[ll]: 



O(ro, d) =< rglro >= exp 
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where Im = (27rpo)~^^^ is the mean spacing between the atoms in the superfluid, and 
the vector d = fq — Fq. The coherent state has the form of eq.(6) below, with / and /' 
replaced by Im, and ji/'e > by the vacuum of the low lying excitations. The above overlap 
integral contains a phase factor, derived from the Berry phase (or Aharonov-Bohm phase 
in this context) of the coherent state. It also contains a gaussian decay factor, reflecting 
the localization of the coherent state. Both factors are characterized by the same length 
scale Im, and are independent of the vortex mass niy. 

In the presence of interactions with the low lying excitations, the total overlap integral 
will change in two ways by the polaron effect[12]: (1) The vortex can induce polarization 
of the excitations, and the overlap between the polarized excitations of one coherent state 
of the vortex and those of a shifted coherent state can contribute to the reduction of the 
total overlap integral. (2) The polarized excitations tend to localize the vortex, squeezing 
the coherent state to a smaller size than /^. These effects will clearly depend on the 
interaction strength, and will also involve the vortex mass. Now a coherent state of the 
vortex centered at tq may be approximated by the following variational wave function 
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where / and /' are two variational parameters, and {t/je > is a wave function of the excitations 
only. With the above ansatz, the total energy of the system is evaluated as 
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where 



H. 



^2. 
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First, the energy is minimized by taking \t(jf, > as the ground state of H^ + Hi, namely 



IV'e >= exp 
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where |0 > is the vacuum of the excitations. The energy of the system then becomes 
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Obviously /' = Im minimizes eq.(lO). The energy is further minimized with respect to / if 



_4 ,-4 , 4m„ [°^ . J{^) ^-k'i' 



r^ = C + ^/ du-^e-'-\ (11) 



'^(^) = E —^ — ^('^^ - ^)- (12) 



where the spectral function JiuS) is defined as 

{k_ 

Having the variational parameters / and /' determined, the overlap integral O(ro, d) is then 
found as 



O(ro,d) =< ro|ro >= exp 






(13) 



for a sufficiently small distance |d|. Here the decay length Id in eq.(13) is 



+ --+/ duj-'-^e-^t.^ . (14) 



1 _ 1 /2 p J(c^) 

The above results have several interesting features. First, the length in the Berry phase 
term is not renormalized by the interactions. In fact, the same result is reached even if we 
assume a general phase factor in the ansatz eq.(6). The result has nicely demonstrated the 
robustness of the Berry phase term against the details of the system. Secondly, by eq.(ll) 
the localization length / is smaller than Z^- The effective mass va^ enters in the equation, 
because it determines the Landau level spacing, which in turn tells how hard it is to mix 
with the higher Landau levels in order to shrink /. Thirdly, the last term of the decay 
length of the overlap integral, eq.(14), comes from the overlap of the polarized excitations. 
Finally, we should point out that when we consider the contribution from the fluctuating 
vector potential in eq.(l) all these features remain unchanged. Eqs. (11-14) are the results 
for the overlap integral from the consideration of the effective theory. 

Now we turn to a completely different way of obtaining the overlap integral, a micro- 
scopic calculation based on Feynman's many-body wave function. We will show that there 
is a complete correspondence between the two approaches. This will enable us to determine 
the vortex mass, the Magnus force, and the coupling to the excitations. If ?/'o(ri...rAr) is 
the ground state many-body wave function of He II, the system with a vortex centered at 
position ro may be described in a first approximation by[10] 

TV 

l'*/'(ro) >= n exp[i6'(rj- - tq) + a{vj - ro)]^o(ri...rAr), (15) 



where ^(r) is the angle of r, and a(r) is a real function of |r|. The most interesting feature 
of the wave function is that it changes phase by 27r whenever an atom moves around the 
vortex center once. In fact, it is by this feature that a vortex state should be defined; 
the above wave function should be regarded as an approximate description of the lowest 
energy state with this feature. The phase factors in eq.(15) introduce a singularity to the 
many-body wave function at the vortex center, and this must be canceled by requiring 
exp[Q;(r)] to vanish at the origin, otherwise the cost in kinetic energy would be too high. 
The particle density in the state, eq.(15), therefore vanishes at yq. At large distances, the 
depletion of particle density due to the vortex vanishes like |r — ro|~^, and correspondingly 
a decays to zero like |r — ro|~^[13]. 

The full calculation of the overlap integral from the many-body wave function is 
difficult, but we may expand lnO(ro, d) in powers of d in the small |d| limit. In this limit 
the two leading terms only involve one- and two-body density distributions in the state 
eq.(15) as will be shown below. Concrete results will then be obtained from a comparison 
with same leading terms in eq.(13). To facilitate the expansion we write 



O(ro,d) =< ij{rQ)\'i{j{ro) >=< exp <j ^[igi{rj -ro,d) + g2{rj -ro,d) 



>, 



(16) 



where <> denotes average over the state of eq.(15), and we have used the notation that 
gi{r,d) = [0(r + d/2)-^(r-d/2)], and that £?2(r, d) = a(r - d/2) + a(r + d/2) - 2a(r). 
Up to second order in d, we may write gi{r,d) = d ■ z x (r — ro)/|r — rop and g2{Y, d) = 
^(d ■ V)^Q;(r). A straightforward cumulant expansion of eq.(16) then yields, to the same 
order in d, that 



lnO(ro,d) = / (frp{r 



+ 



- I (frp{r) 



.d ■ z X (r — fq) 

' i 12 

d ■ z X (r — To 



+ / d^rp{r) 

1 2 



l(d.V)Mr-ro) 



d rd r p(r, r 



r-ror 
, d ■ z X (r — To) d ■ z X {r' — tq] 



(17) 



d'^Y p{v] 



d- z X [r — fq) 



r-ro| 

H 2 



I / 19 

r — To 



Ir-rol 



where p(r) =< ^ . d{r — Tj) > and p(r, r') =< J2i=ii ^(^ ~ rj)5(r' — r_,) > are the one- and 
two-body density distributions in the state eq.(15). 

The first order contribution to lnO(ro, d) in eq.(17) is purely imaginary, which can be 
evaluated as iirpoz ■ (d x fq) if we replace p(r) by po, assuming that our system (including 
the vortex center) is confined within a disc centered at the origin of r. Here po is the 2-d 
superfiuid number density. The correction due to p{r) — pQ = pi{r) is zero in the infinite 
size limit, because of the rotational symmetry in pi{r) about ro and the fact that the 
density depletion decays sufficiently fast at large distances from tq. This first order term 
is the Berry phase associated with the Magnus force discussed in Ref.[14]. 

The second order contribution to lnO(ro, d) in eq.(17) is purely real, and it must also 
be negative as required by the fact that |0(ro,d)| < 1 for nonzero d. We may therefore 
write 



lnO(ro, d) = iwpoz ■ (d x fq) — —^ + higher order terms, (18) 

where we have put the second order term as independent of the direction of d because of 
the isotropy of the system about the vortex center. The second order coefficient has been 
parameterized by Id, which has the dimension of a length, and represents the same decay 
length as in eq.(13). 

We now examine closely the second order terms in eq.(17), and show that their con- 
tribution to /^ is finite. The term containing a{r — ro) converges because the double 
derivative of a decays as an inverse cubic function at large distances from the vortex cen- 
ter while p{r) approaches a constant. At short distances, a may diverge like a logarithm, 
but p{r) vanishes linearly, causing no trouble to the convergence of the integral. Therefore 
we shall no longer consider this term. In the presence of particle correlation, the form 
of p{r, r') is unknown for the state containing a vortex, except at large distances away 
from ro, where it reduces to po(r — r'), the distribution in the absence of the vortex. We 
may, however, replace the distributions by their asymptotic forms in eq.(17) in order to 
examine the long distance contributions to these terms, because it is only from there that 
a divergence may ever be possible. Then, the last three terms of eq.(17) (c./. eq.(18)) yield 



where '...' stands for the correction due to short distance contributions, S'o(k) is the static 
structure factor in the absence of the vortex, and F(k) = i2Tve'^^'^"z-dx k//c^ is the Fourier 
transform of d • i x (r — ro)/|r — rop. It is known[13] that 5'o(k) = hk/2Mc for small /c, 
where M is the mass of a helium atom, and c is the sound velocity. The integral in eq.(19) 
therefore converges, meaning that the second order expansion in eq.(17) exists in realistic 
situations. 

Before we proceed further we would like to comment on the validity of the above 
discussions. We have ignored multi-particle correlations induced by the vortex in the 
original wave function eq.(15) and in the evaluation of the expression eq.(17). We assume 
that the induced correlations decay sufficiently fast away from the vortex center, such 
that they do not affect the convergence properties at large distances. The situation at 
short distances is very complicated[10], and the short distance contribution can be quite 
substantial to the reduction of the overlap function. We expect, however, that the system 
should behave smoothly at short distances, so that no divergence can be induced from 
there. Our later arguments will only be based on the conclusion drawn above that the 
decay length /^ is finite. 

The decay length l^ strongly depends on the interaction between the atoms in the 
superfluid. As the interatomic interaction becomes weak, the sound velocity decreases, 
which makes 5*0 (k) large and therefore Id small from eq.(19). In the extreme case of no 
interatomic interaction, /^ becomes zero. This is just what one should expect from a direct 
calculation of eq.(17) in the free boson limit, in which case p(r, r') = p{r)p{r'){N — 1)/N. 

With the overlap integral evaluated both from the effective theory, eq.(13), and Feyn- 
man's many-body wave function, eq.(18), we now would like to see how the parameters of 
the effective theory should be constrained. Firstly, the parameter po in the Magnus force 
of the effective theory is the same as the 2-d superfluid density from the comparison of the 
results for the Berry phase term of the overlap integral. Second, in order to be consistent 
with the result of finite Id from Feynman's many-body wave function, eq.(19), the sum in 
eq.(14) has to be convergent, implying that the spectral function J{uj) must vanish faster 
than u) at low frequencies and |M(/c)| must be less singular than k~^ at small /c's. A 



comparison of eqs.(19) and (14) suggests that in the low frequency Umit 

Ji^) - ^y. (20) 

In the language of quantum theory of dissipation[15], this kind of coupling is of the so- 
called superohmic type. In a recent study of vortex tunneling in Ref.[16], a general heat 
bath is considered. It is found there that a superohmic coupling to the heat bath has a 
negligible effect on the tunneling process at low enough temperatures. 

As for the mass of the vortex, our result of finite decay length implies that the mass of 
the vortex cannot be infinite, otherwise the localization length of the vortex would shrink 
to zero according to eq.(ll) and the decay length of the overlap function would become 
zero according to eq.(14). Therefore, our result is consistent with that of Refs.[6,7], which 
suggest that rriv is zero or finite, and is in apparent disagreement with that of Ref. [8] 

The vortex mass that we originally introduced in eq.(l) may have already included 
the effect of renormalization by the polarization of all but the low lying excitations of the 
superfiuid. There is still a possibility that it may be renormalized to infinity if the polar- 
ization of the low lying excitations is included. Indeed, if we neglect the Magnus force, a 
straightforward perturbative calculation[12] shows that the mass renormalization becomes 
logarithmically divergent if the coupling spectrum J{uj) goes as u)"^ at low frequencies. This 
is essentially the result in Ref. [8]. The divergence of the mass renormalization becomes 
severer if J{oj) would vanish slower than uj"^. However, this will not be the case for a 
dynamical process with a time scale, as the following arguments show. 

The situation in the presence of the Magnus force is quite different. One can no 
longer set up a momentum eigenstate and extract an effective mass of the vortex from the 
energy dependence on the momentum. A more natural approach is to relate the effective 
mass to higher Landau levels of cyclotron frequency Uc = hpo/niy. Interaction with low 
lying excitations may shift and broaden the higher Landau levels, but these effects are 
not divergent in a perturbative calculation using eq.(2) if the coupling is superohmic. 
Therefore, if the higher Landau levels are well defined before turning on the coupling to 
the low lying excitations, we can conclude that further inclusion of such coupling has little 
effects on the higher Landau levels and thus the effective mass of the vortex. To observe a 

9 



higher Landau level experimentally, one may trap ions in vortices produced in a rotating 
film of superfluid, and excite the vortices by electrical coupling to the ions. [17] 

Finally, we would like to make some remarks about the generality of our results. As 
long as the Feynman many-body wave function description of the vortex state is valid, 
everything else just follows from standard many-body physics such as the form of 5*0 (k) at 
small k. As long as 5*0 (k) vanishes with some positive power of k, orthogonality catastrophe 
in the overlap integral will not occur. Therefore, our results may also be applicable to 
vortex structures in superconducting films and wire networks, Josephson junction arrays, 
and quantum spin systems. 
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